X-ray Raman Scattering from Water Near the Critical Point: 
Comparison of an Isotherm and Isochore 
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X-ray Raman spectra of liquid, sub- and super- critical water at the oxygen K-edge were measured, 
at densities p — 1.02 - 0.16 gem -3 . Measurements were made along both an isotherm and an isochore 
passing near the critical point. As density is reduced there is a general tendency of the spectra to 
increasingly resemble that of the vapor phase, with, first, a well separated low-energy peak, and, 
eventually, at densities below the critical density, peaks appearing at higher energies corresponding 
to molecular transitions. The critical point itself is distinguished by a local maximum in the contrast 
between some of the spectroscopic features. The results are compared to computed X-ray absorption 
spectra of supercritical water. 

PACS numbers: 78.70.Ck, 33.20.Rm, 61.25.Em, 31.15.A-, 33.15.xv 
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Heating a liquid under pressure allows one, eventually, 
to pass into the super-critical region where there is a 
continuous change in density between liquid and vapor 
phases. The neighborhood of the last point where there 
is a discontinuous change between the two phases, the 
critical point (see Figdjleft)), is interesting for a variety 
of reasons. Fundamentally, the critical point is where 
strong density fluctuations occur on essentially all length 
scales, leading to unique physical properties. Practically, 
the properties of liquid solvents in the supercritical region 
facilitate different types of reactions, often without the 
(sometimes environmentally harmful) catalysts that can 
be needed in other conditions [H-13]. 

Water in the critical and supercritical region (critical 
point, T c = 647.096 K, P c = 22.064 MPa, p c = 0.322 
gem -3 [5J]) has attracted significant attention both for 
fundamental and practical reasons. On the fundamental 
side, there is great interest in the nature of the bonding 
in this region: how do the hydrogen bonds, which are 
responsible for many of the interesting features of liq- 
uid water in ambient conditions, survive into the super- 
critical region? There have been several studies focusing 
on this topic (e.g. experimentally @-18| and theoreti- 
cally [H-tlBi)- 

However, there remains significant debate 
and uncertainty, in part because experiments on water 
in this region (where, for example, its reactivity is suffi- 
ciently high to etch typical stainless steel) are difficult, 
and in part because the interpretation of spectroscopic 
experiments is severely complicated by the density fluc- 
tuations. 

In the present paper we apply X-ray Raman Scattering 
(XRS) to investigate the behavior of water in the critical 
region and compare our results to ab-initio calculations. 
Soft X-ray absorption spectroscopy (XAS) [13, 28 1 and 



XRS [29( have revealed the details of hydrogen bonding 
(HB) environments in water by assigning the observed 
spectral features near the oxygen K-edge 17, 3p|- 35l| . 
(Non-resonant) XRS is alternative to XAS. Under the 
condition qro <C 1 (ro: radius of core state), the dipole 
contribution is dominant and has the same matrix ele- 
ment as absorption. In this case, the spectral shape of 
XRS is proportional to XAS cross section. The higher 
energy of XRS experiments (6 - 12 keV), as compared to 
XAS work (~ 540 eV), means that it is much easier to 
penetrate into complex sample environments, as is very 
important for investigating supercritical water (SCW). 

The experiments were performed on the Taiwan beam- 
line BL12XU (3(| at SPring-8 in Japan. Incident ra- 
diation was monochromatized by a Si(400) reflection 
(—,+,+,—) 4-bounce high-resolution monochromator 
and focused onto the sample to a spot size of 80 (V) x 
120 (H) /im 2 . The scattered radiation is analyzed using a 
2-m radius spherically bent Si(555) analyzer crystal in a 
Rowland circle geometry. The analyzer angle of was fixed 
at 88.5°, corresponding to an analyzer energy of 9888.8 
cV. The incident photon energy was scanned from 10421 
to 10437 eV. The total energy resolution was 260 meV 
(FWHM of the quasi-elastic line of the sample) measured 
at the analyzer energy. The energy scale was calibrated 
using scans of the Tantalum L3 (9881.1 eV) and Rhenium 
L 3 (10535.3 eV) absorption edges. 

The sample cell was custom designed for inelastic X- 
ray scattering experiments, with an inner chamber made 
of Hastelloy-X alloy with diamond windows to avoid reac- 
tion with the SCW. For this work, we chose a 3 mm sam- 
ple length, corresponding to the 1/e absorption length for 
water of density 0.6 gem -3 . The incident and outgoing 
angular acceptances were 22°. This cell was then placed 
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inside of a vacuum chamber (to aid thermal control and 
reduce air-scatter) with polyimide windows. 

The sample was ultra-pure water. Its thermodynamic 
conditions were maintained by heating the inner sample 
cell and providing pressure using an external hand press. 
As the sample density is extremely sensitive to temper- 
ature and pressure especially near the critical point, the 
pressure was continuously monitored throughout the ex- 
periment and the temperature was controlled using an 
Inconel-covered thermocouple in direct contact with the 
water in the cell (just to one side of the X-ray path). 
Both pressure and temperature measurement systems 
were carefully calibrated, with the temperature gauge ex- 
pected to be accurate to 0.4 % and the pressure to ± 0.1 
MPa and 0.1 %. During scans, measured temperatures 
were maintained to ± 0.1 K and pressures to ± 0.05 MPa. 

Measurements were made at an 18° scattering angle, 
or q ~ 15.7 nm _1 . This satisfies the dipole scattering ap- 
proximation with qr\ ~ 0.1 <C 1 (ri = clq/Z is the radial 
extent of the oxygen Is wave function, where, ao is Bohr 
radius and Z is effective nuclear charge for the orbital). 
The acceptance of the analyzer was 2.5 nrn^ 1 . Scans 
over the full energy range, 532 - 548 eV energy transfer 
were typically 90 minutes long (including a check of the 
elastic peak) and depending on the conditions, e.g. sam- 
ple density, typically 10 to 20 scans were measured at 
each set of thermodynamic conditions, with, sometimes, 
a shorter range chosen to efficiently use beam time. For 
each scan, the signal was normalized by the incident in- 
tensity and backgrounds were subtracted assuming an 
exponential energy dependence, i.e. primarily, due to 
the Compton tail of the diamond windows. The spectra 
were then normalized to unity integral over 532 to 548 
eV energy transfer, with a smooth continuation used if 
the measured spectra did not cover the full range. 

The measured conditions arc indicated on the phase 
diagram in Figfl] with the precise conditions given in 
Table Q] FigJSjA) and (B) present the measured spec- 
tra. The XRS/XAS spectra can be divided into three 
regions: (I) pre-edge (533-535 eV), (II) main-edge (535- 
538 eV), and (III) post-edge (> 538 eV). Fig^A) shows 
the isothermal response at the critical temperature and 
FiglJjB) the response along isochore at the critical den- 
sity, both passing close to the critical point (spectrum 
(d)). The clearest trend is visible along the isothermal 
line, showing a gradual progression from the relatively 
blurred out spectra at high density (b), similar to bulk 
water, to something closer to the gas-phase spectra at 
low density (f). There is also a shift, by about 0.5 eV, 
to lower energy of the pre-edge when the temperature is 
increased from 293 K to 653 K (spectrum (a) to (b)). 
The lowest energy peak corresponds, in the gas phase, 
to the Is — > 4ai transition of the H2O molecule, which 
has an excited state wave function spread out over both 
hydrogen atoms and separate portion about the oxygen 
atom. This peak becomes more distinct with decreasing 
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FIG. 1. Left: (T-P) phase diagram and isochores of H2O 
5]. The solid line corresponds to the liquid-gas coexistence 
curve. Dashed lines indicate the expanded process and solid 
dots indicate the measured conditions. Right: (p-P) phase di- 
agram and isotherm of H2O. The dashed line indicates critical 
isochore. 



TABLE I. Thermodynamic conditions and uncertainties at which 
the XRS spectra of expanded water were measured. The spectra 
are recorded along one near-critical isotherm T ~ 1.01 T c (b, c, d, 
e, f) and along one critical isochore p = p c (g,h,d). Labels indicate 
in Fig. [1] 





T (K) P (MPa) 


p (1) (gcm 3 ) 


P+ 


P- 


p (2) (gcm 3 ) 


a 
b 
c 
fl 
e 
f 

g 
h 


293.15(1.6) 54.6(0.15) 
653.15(1.6) 54.6(0.15) 
653.15(1.6) 29.6(0.15) 
653.15(1.6) 23.3(0.15) 
653.15(1.6) 22.8(0.15) 
653.15(1.6) 22.0(0.20) 
759.25(2.0) 54.0(0.20) 
708.15(1.8) 39.1(0.20) 


1.02(0.00) 
0.64(0.01) 
0.53(0.01) 
0.29(0.09) 
0.21(0.03) 
0.17(0.01) 
0.32(0.01) 
0.32(0.01) 


1.02 1.02 
0.64 0.65 
0.52 0.54 
0.21 0.34 
0.18 0.25 
0.16 0.17 
0.32 0.33 
0.31 0.33 


1.02(0.00) 

0.64(0.00) 

0.53(0.01) 

0.24(±°;Jt) 

0.20(^° 6 2 ) 

0.32(0.01) 
0.32(0.01) 



a T: nominal temperature, P: nominal pressure, pM: center of 
the error value at T, p_|.: upper limit of the error at T+AT, 
p_: lower limit of the error at T-AT, p( 2 ': nominal value 



density or pressure. Meanwhile there are also changes in 
the high-energy part of the spectrum, but these are hard 
to quantify, especially in the absence of a rather good 
model. The trends in the isochore, varying from 1.01 to 
1.19 T c , are not so strong. 

The spectrum at the critical point, (d), is distinguished 
from the others by having a relatively strong contrast be- 
tween the low-energy (presumably Is — ¥ 4ai) peak and 
the high-energy part of the spectrum. This is brought 
out in more detail in Figj3l where the difference between 
the maximum and minimum intensities is plotted. The 
contrast in spectrum (d) is stronger than that in spec- 
trum (c) and (e) (at T c but, respectively, at higher and 
lower pressure) and (h) on the high-temperature side of 
the critical isochore. This contrast is due to an increase in 
height of the low-energy peak relative to the high-energy 
edge. If one considers the response in the critical region 
to be the sum over the spectral response of clusters of 
various sizes, improved contrast in the conditions where 
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FIG. 2. (a)-(h): Area normalized XRS spectra in expanded 
water; near K-edge spectra of the oxygen obtained with a to- 
tal energy resolution of 260 meV (FWHM) and 0.2 eV step. 
Open circles with statistical error bars are measured data and 
solid lines are smoothed curve. The thermodynamic condi- 
tions are indicated in Table [I] (Gas Phase): Dilute gas phase 
spectrum (fluorescence-yield soft-XAS from Ref.[37]]) convo- 
luted by resolution function of this work. (A) along with 
isotherm T = 1.01 T c (B) along with critical isochore p — p c . 
(C) Computed XAS spectra of SCW for several densities at 
T = 653 K. See text and supplemental material. 
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FIG. 3. Density dependent contrast of pre-edge with maxi- 
mum near the critical point. 



the largest distribution of cluster sizes is expected is sur- 
prising. Thus these measurements might indicate that we 
are seeing effects from electronic motions on time scales 
comparable to the XRS. 

Our data is noticeably different than a recent publi- 
cation [13] of data in nominally similar thermodynamic 
conditions. This difference might be the result of a large 
non-dipole contribution in Ref.[17| (they had qr\ ~ 0.45) 
or other experimental differences (they had somewhat 
worse, ~1 eV energy resolution and oxygen- containing 
sapphire windows). If we apply their scaling argument, 
the size of the peak at 534 eV suggests about 73(10) % 
of the molecules are nearly gas-like in spectrum (c) (th e 
point that was closest to the conditions stated in Ref.[l7j) 
and about 85(10) % in spectrum (d), a much-larger gas- 
like percentage than the 35(20) % of Ref.[I3|. How- 
ever, such a scaling argument probably over-simplifies 
the problem and detailed simulations are really required 
to interpret the results. 

Simulations were performed using Car-Parrincllo 
molecular dynamics (CPMD) [38| in the microcanonical 
ensemble with density functional theory. Initial struc- 
tures were obtained from the DL_POLY (39[ simulations 
for densities 0.73, 0.64, 0.55, 0.38 and 0.35 gem" 3 at 653 
K using the TIP2P water model. Exchange-correlation 
function [40I ] and norm-conserving pseudopotentials [ill ] 
were used. Oxygen if -edge XAS spectra were computed 
using the transition state potential approach of Slater 
[42| for the geometries obtained in the molecular dynam- 
ics runs. Details of the computation are presented in the 
supplementary materials. 

The calculated XAS along isotherm T = 653 K are 
shown in Figl^C). A comparison of the calculated XAS 
of SCW and vapor-phase from experiment shows the 
presence of the pre-edge peak in all of the spectra at 
a similar excitation energy (534 eV). The calculated pre- 
edge peak height increases with the decrease in density 
consistent with the trend of experimental results. The 
main-edge region narrows with the decrease in density 
in both experiment and calculation, though this is not 
clearly consistent especially near the critical density. 

Examination of the final state wavefunction (Figf4]) 
indicates that the origin of pre-edge peak (I) is excita- 
tions to final states that are a combination antibonding 
OH and Rydberg orbitals. Examination of the orbital 
plots shows that these states are confined to the wa- 
ter molecule with the tail of the wavefunction extend- 
ing only to nearest neighbor molecules. Thus, changes 
in the HB at the excited water molecule will have little 
effect on the shape and energies of these highly local- 
ized orbitals. This explains why the pre-edge peak posi- 
tion and profile are largely unaffected by the surrounding 
environment, showing no significant modification with a 
change in the density of SCW. Examination of orbital 
plots from the calculated XAS indicate that transition 
in the main-edge region (II) are to final states with an- 
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FIG. 4. Molecular orbital plots for the final states with re- 
spect to the three excitation regions (I: pre-edge, 533-535 eV; 
II: main-edge, 536-538 eV; III: post-edge, > 538 eV) . The two 
thermodynamic conditions, p = 0.64 gem -3 (= 2.0 p c ) and p 
— 0.35 gem -3 (= 1.09 p c ), are selected. H2O molecules are 
shown as v-shaped sticks with ball. The atoms color code is 
red for O and light blue for H. The molecular orbitals are plot- 
ted in positive (pink) and negative (purple) valued wavefunc- 
tions using isosurface values of (region I) 0.05 eA -3 ; (region 
II) 0.02 eA -3 ; (region III) 0.015 eA~ 3 . 



tibonding OH orbitals and a larger mixture of Rydberg 
states than found in the lower-lying excitations that are 
responsible for the pre-edge peak (I). At lower densities, 
where less interactions with surrounding water molecules 
are expected, there is a lower mixing in the diffuse Ryd- 
berg states, resulting in an increase of the energies and 
a narrowing of the distribution of these states. This is 
the reason for the narrower main-edge peak in the 0.35 
gcm~ 3 spectrum compared to 0.64 gem -3 . An examina- 
tion of the post-edge region (III) orbital plots for SCW 
indicates that excitations in this region are to diffuse, 
continuum orbitals. A decrease in density shows some 
localization of the wavefunctions of these states between 
the molecules. Both measurements and calculation show 
sensitive changes to this region promising additional in- 
formation if calculations are improved. 

In conclusion, measurements of the oxygen K-edge x- 
ray Raman scattering show marked changes as thermo- 
dynamic conditions are tuned in the neighborhood of 
the critical point. There is, broadly speaking, reason- 
able qualitative agreement with XAS calculations that 
allows us to interpret these results in terms of orbital 
structures. However, the contrast in the measured data 
is better than in the calculation, especially near the ex- 
act critical conditions, which hints that we may be seeing 
an effect of fast electronic dynamics in the XRS spectra. 
The complex high-energy structure of the spectra and 
calculations also suggest significantly more information 
might be obtained if the calculations can be modified to 
bring them into closer agreement with the data. 
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